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ABSTRACT 

We present the optical and X-ray properties of four clusters recently discovered by the South Pole 
Telescope (SPT) using the Sunyaev-Zel'dovich effect (SZE). The four clusters are located in one of 
the common survey areas of the southern sky that is also being targeted by the Atacama Cosmology 
Telescope (ACT) and imaged by the CTIO Blanco 4-m telescope. Based on publicly available griz 
optical images and XMM-Newton and ROSAT X-ray observations we analyse the physical properties 
of these clusters and obtain photometric redshifts, luminosities, richness and mass estimates. Each 
cluster contains a central elliptical whose luminosity is consistent with SDSS cluster studies. Our 
mass estimates are well above the nominal detection limit of SPT and ACT; the new SZE clusters are 
very likely massive systems with M^5x 10 14 M . 

Subject headings: cosmic microwave background — cosmology: observations — galaxies: distances 
and redshifts — galaxies: clusters: general — large-scale structure of universe 



1. INTRODUCTION 

Galaxy clusters serve as important cosmological probes 
as their formation and evolution rate depend on cosmo- 
logical parameters and the kinematics of the dark matter. 
They are the observational counterparts of dark matter 
halos and their abundances and masses as a function of 
redshift are quite sensitive to the growth of structure in 
the Uni verse offering a potenti ally powerful probe of dark 
energy ( Carlstrom et al.|2002[ ). Galaxy clusters also har- 
bor a significant fraction of the visible baryons in the 
Universe. Hydrogen gas, heated as it falls into the dark- 
matter potential well of a cluster, forms an intraclus- 
ter medium too hot to be bound by individual galaxies. 
The hot intracluster gas leaves an imprint on the Cosmic 
Microwave Background Radiatio n (CMBR) though the 
Suny aev Zel'dovich effect (SZE) (Sunyaev & Zel'dovich 



1972 



terec 



in which CMB photons are inverse-Compton scat- 
by the hot intracluster gas. It is this SZE sig- 
nal that two new ground-based mm-ban d telescopes, the 



Atacama Cosmolo gy Telescope (ACT) (Kosowsky 



2006 



Fowler et al | |2007| ) and the South Pole telescope (SPT, 
'( Ruhl et al.|2004p , have been designed to detect through 
its frequency-dependence: these experiments will mea- 
sure intensity shifts shifts of the CMB radiation corre- 
sponding to a decrement below and an increment above 
the "null" frequency around 220 GHz. Both telescopes 
are now actively acquiring data over large areas of the 
southern sky that will ultimately cover several hundreds 
to thousands of square degrees. 

Recent results f r om the SPT collaboration 
(|Staniszewski et ah] [2008]) have yielded the first 



blind detection of SZE clusters in one of the common 
southern survey regions centered near right ascension 
05 lu '30 m and declination -53°. This area has also been 
scanned by ACT and has been optically imaged in the 
optical (griz) with the CTIO Blanco 4-m telescope as 
part of the Blanco Cosmology Surve y (BCS) (Menanteau 
et al |2008||Staniszewski et al.|2008| ). Prompted by these 
results, here we present a detailed study of the physical 
properties of these first four SZE-detected clusters based 



on publicly available optical imaging and archival X-ray 
data. Throughout this paper we assume a flat cosmology 
with H a — lOOh km s _1 Mpc -1 , h — 0.7 and matter 
density Am = 0.3. 

2. SZE OPTICAL COUNTERPARTS 

The positions of the four SZE-selected clusters re- 
ported by the SPT team are contained within the region 
surveyed by the BCS in 2005 and 2006. These obser- 
vations have been publically available for about a year 
or more now. We have processed and analyzed them 
using an independent software pipeline develop ed by us 
at Rutgers Unive rsity. We refer the reader to |Menan-| 



teau et al. ( 2008 1 where we provide a full description of 
the observing strategy, pipeline and data products and 
some results on new massive galaxy clusters discovered 
in the BCS survey field centered near 23hr. Here we use 
the identical software pipeline to generate data products 
from which we characterize the optical properties of the 
four SZE clusters. First, we identify the optical coun- 
terparts, which appear as clearly visible overdensities of 
early-type galaxies near the sky locations of the SZE 



decrements (see Table 1 from Staniszewski et al. 2008 1 



Roughly centered within each overdensity is a bright el- 
liptical, which we take to be the Brightest Cluster Galaxy 
(BCG) . In Figure [T] we show the gri color images of the 
four systems, centered on the location of the BCG. Each 
system shows a dominant population of early-type galax- 
ies with very similar colors that can easily be identified 
by visual inspection. In all cases the offset between the 
BCG and the SZE decrement, shown as red ellipses and 
green crosses respectively, is less than ~ 30 arcsec. 

2.1. Photometric Redshift Determination 

We determine photometric redshifts fr om th e four- 
band optical images using BPZjBem'tcz 200C| follow- 



ing the same procedure as in |Menanteau et al"| ( |2008[ ) 



To avoid contamination by non-cluster members we use 
the BCG to estimate the photometric redshift of each 
system. The BCG is, by definition, indisputably part 
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3PT-CL 0528-5^1; z=0.7l> 
• 




S0520 ( |Abell et al.|1989) which is also the X -ray cluster 
RXC J0516.6— 5430 ( jBohringer et al I [2004b Its pub- 
lished spectroscopic redshift of z = 0.294 dGuzzo et al. 



1999 1 is in excellent agreement with our photometric red 
shift for that BCG, z = 0.27 ± 0.02, but s omewhat dif- 
ferent from the value of z — 0.35 quoted by |Staniszewski| 



Fig. 1. — Compos ite pri color image of the four SPT SZE-selected clusters from the data. Green crosses show the position of the 
SZE decrement from [Staniszew ski et a,l. | ( |2008| ) for each system. The red ellipses denote the position of the BCG and show as well the 
SExtractor Kron radius. For cluster yFT-CL U517— 5430 (Abell 0520, main panel) we show the XMM-Newton X-ray isointensity contour 
levels (starting near the center at a value of 1 X 10~ 4 cts s —1 pixel -1 and decreasing going outward by factors of 2.5). 

of the cluster, resides in a quasi-central location and al- 
ways provides the strongest signal and the best colors 
since such galaxies are generally several times brighter 
than L* at any given redshift. Moreover, the luminosities 
and particularly the colors of BCGs are well constrained 
as they are very bright elliptical galaxies dominated by 
old metal-rich stellar populations for which multi-color 
spectral energy distribution (SED) fitting techniques are 
well suited yielding precise and reliable results. In the 
bottom panel of Figure [2] we show the output BPZ pho- 
tometric redshift probability density function (pdf) for 
each BCG. The differences between the BPZ Bayesian 
and Maximum Likelihood (ML) predicted redshifts are 
small (8z ~ 0.01) in all cases and therefore we will em- 
ploy ML photometric redshifts hereafter as the addition 
of a prior is not justified when, as is the case here, the 
BCG SEDs are fitted unambiguously. We note that oc- 
casionally BCGs host radio- loud active galactic nuclei 
(Best et al. 20071. Although we have no indication of 
this, if any of our BCGs do harbor AGN then we would 
expect some blueing of the colors causing our photo-z's 
to be slightly over-estimated. 

Only one of the clusters, SPT-CL 0517-5430, has been 
previously reported; it is the optically-rich cluster Abell 



et al. ( |2008[ ) for the cluster which came from the peak of 
their photo-z distribution of red-sequence galaxies using 
the same BCS imaging. In Table I] we present the clus- 
ters' optical positions based on the BCG locations and 
photometric redshifts with l-a uncertainties. 

2.2. BCG luminosities 

The range in brightness of BCGs has been observa- 
tionally established from galaxy clusters in the SPSS up 
to redshifts of 0.36 (see, for example, Loh & Strauss 



2006| |Koester et al.||2007[ ). This can be used to deter 
mine whether the BCGs we identify in the SZE clusters 
have luminosities consistent with SDSS BCGs. In Fig- 
ure [2] (2nd panel from top) we show the total r-band ob- 
served magnitude of each cluster's BCG as a function of 
redshift compared to a parametrization of the observed 
r-band magnitudes of SPSS B CGs (dashed line) (Fig- 
ure 7 in Loh & Strauss 2006). The curve corresponds 
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Fig. 2. — The absolute (upper panel) and observed (2nd panel) 
r-band magnitude of each BCG as a function of their photomet- 
ric redshift. The dashed line represent s the modeled observe d r- 
band magnitude of BCGs according to |Loh &: S trauss ( 2006JI . In 
the third panel we show the clusters sort X-ray tlux as a function 
of redshift. The solid and dashed curves represent the expected 
X-ray fluxes for clusters with mass of M200 = 1 x 10 15 Mq and 
3 X 10 14 Mq respectively. The bottom panel shows the photo-z 
probability distribution function for each BCG and with a cross 
at the location in the curve for SPT-CL 0517— 5430's spectroscopic 
redshift. Error bars in the figure represent Her deviations drawn 
from the probability functions. 

to t he simple prescription M* — 1.5 where M* is taken 



from Blanton et al. ( 2003 1 and the model includes passive 
evolution with redshift. We have extrapolated the SDSS 
BCG luminosity model beyond z — 0.36 as there is over- 



sively evolved since z ~ 1 1 


Blakeslee et al. 2003; Lidman 


let al. 12008 Mei et al. 112008 


I . We see from the figure that 



model curve (i.e., are intrinsically brighter) and we take 
this as strong evidence in support of them being bona 
fide BCGs. It is worth noting that SPT-CL 0509-5342 
at z = 0.36 is the least luminous of the four, just slightly 
brighter than the luminosity predicted by the model. As 
we will see in the next section this cluster also has a 



lower mass and X-ray luminosity compared to the other 
systems. In Figure 2] we also show the range of absolute 
magnitudes of the BCGs as a function of redshift. Here 
we again see that SPT-CL 0509—5342 is the least lumi- 
nous object, ~ 1 magnitude fainter than the other BCGs 
which all have about the same luminosity. 

3. CLUSTER MASS ESTIMATIONS 

The new SZE experiments aim to deliver mass-selected 
samples of galaxy clusters out to large redshifts which can 
be used to probe the growth of structure in the Universe 
and to study cluster phys ics. Current N-body plus hy- 
drodynamics simulations ( Motl et al.||2005 Nagai 2006 



Sehgal et al.| 2007[ ) predict robust scaling relations be 
tween the SZE signal (i.e., the y-distortion due to inverse 
Compton scattering) and cluster mass. However, a fun- 
damental step before using SZE clusters as cosmological 
probes is to observationally establish the relationship be- 
tween y and mass using mass estimates independent of 
the SZE measurements. SZE-sclected clusters will almost 
surely be the focus of intense efforts to obtain mass esti- 
mates using X-ray images and spectra, weak and strong 
lensing and velocity dispersions. In the following, we 
present the first independent mass estimates for these 
new SZE-sclected clusters from archival optical and X- 
ray data. 

3.1. Optical Cluster Masses 

Here we use optically observed parameters of galaxies 
to predict the cluster mass from weak lensing ric hness- 
mass scaling relations established f rom the SDSS ( John 



ston et al.||2007| |Reyes et al.||2008 ). We a pply the state 
of-the-art mass tracer (Reyes et al. 2008) built from the 



SDSS maxBCG cluster sample ( |Koester et al.||2007| ) of 
~ 13,000 optically-selected clusters. Here we briefly de- 
scribe our procedure and refe r the reader to our earlier 
study ( Menanteau et alT]|2008 ) for complete details. The 

l.oquantities, -/V|q , L200, and £bcGj are the observational 
values needed as input to the cluster mass scaling rela- 
tion. The cluster richness, N^qq, is the number of E/S0 
galaxies within i?2oo with colors and luminosities that 
satisfy specific conditions for membership. i?200 is the 
estimated cluster size defined as the radius where the 
cluster galaxy density is 200il J „ 1 times the mean space 
density of galaxies in the present Universe. Similarly, 
L200 is the total rest-frame integrated r-band luminos- 
ity of all member galaxies included in N^qq and Lbcg 
is the rest-f rame r-band luminosity of the BCG. |Reyes| 



et al. ([2008]) provide power-law functions for both the 
luminosity-mass and richness-mass relations (see section 
5.2.1 in their paper), although here we only use the re- 
lation based on luminosity, which should be more robust 
than that based on richness. In Table [T] we present the 
cluster mass estimates and luminosities where M(L2oo) is 
the mass observational equivalent of A^oop 1 ■ For the two 
higher redshift clusters, SPT-CL 0528-5300 at z = 0.70 
and SPT-CL 0547-5345 at z = 0.88 we quote lower-limit 
estimates as some fraction of low luminosity galaxies fall 
below our magnitude- limit (i ~ 22.5). If we make the as- 
sumption that the two higher redshift clusters have lumi- 

M20OP is the halo mass enclosed within a radius of spherical 
volume within which the mean density is 200 times the critical 
density. 
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nosity functions that are similar to SPT-CL 0517—5430 
at z = 0.29, then we estimate that we arc missing 15% 
and 58% of the luminosity for SPT-CL 0528-5300 and 
SPT-CL 0547—5345 respectively. The mass estimates in 
Table [I] include these corrections. 



3.2. X-ray Properties and Masses 

The nearby cluster SPT-CL 0517-5430 has accurate 
measurements of its X-ray luminosity (Lboi = 9.2 ± 1.2 x 
10 44 ergs s _1 ) and temperature ( kT = 7.5 ± 0.3 keV ) 



from XMM-Newton observations ( |Zhang et al.| [2006 1 . 
There is also a hydrostatic mass estimate of -M500 = 
6.4 ± 2.1 x 10 14 Af Q (note this is at an overdensity of 
500 times the critical density). The cluster's X-ray mor- 
phology is strongly elongated: it has an axial ratio of 
1.4-1.7 with its major axis aligned ^14° from north (see 
Fig-Q. 

We utilize the ROSAT All Sky Survey (RASS) to ob- 
tain X-r ay information on the oth er three clusters. As 
noted by Staniszewski et al. (20081, faint RASS sources 
lie close to the locations of SPT-CL 0509-5342 and SPT- 
CL 0547—5300. For purposes of determining X-ray fluxes 
and luminosities we will assume that these RASS sources 
are indeed the X-ray counterparts to the SZ clusters. 
Starting with the raw X-ray photon event lists and ex- 
posure maps from the MPE ROSAT Archive 2 we deter- 
mined count rates for the band covering PI channels 52 
to 201. Extraction regions were optimized to obtain all 
detected source photons (a radius of 3' was sufficient) 
and, for the background estimation, to reduce statistical 
fluctuations (using an annular region between 5' and 25'). 
The X-ray source position was used for the two SZ clus- 
ters with apparent counterparts; for SPT-CL 0528—5300, 
the quoted SZ position was used to derive upper limits. 

The background-subtracted count rates (SPT- 
CL 0509-5342: 0.025 cts/s; SPT-CL 0528-5300: 
< 0.015 cts/s; SPT-CL 0547-5345: 0.007 cts/s) were 
converted to fluxes assuming a hard thermal spectrum 
(kT ~ 5 keV) and accountin g for each cluster's Galacti c 
absorbing column density (Dickey & Lockman 19901; 
values fell in the range 6.1 x 10 zu atoms cm"' to 7.8 x10 
atoms cm~ 2 for all sources. A fc-correction was also 
applied to convert fluxes to the standard 0.1-2.4 keV 
ROSAT band in the rest frame of each cluster. These 

(2nd panel from bottom). 
0528-5300 is at 2-cr and 



fluxes are plotted in Fig. [2 
The upper limit for SPT-C 
the soft 0.1-2.4 keV band flux for SPT-CL 0517-5430 
was converted from the bolometric luminosity quoted by 
Zhang et al. (2006). We denote the region (at the high 



Source Catalog from the RASS. Also shown are curves 
of the expected X-ray fluxes for clusters with masses of 
A/200 = 1 x 1O 15 M (solid) and 3 x 10 14 Af Q (dotted). 
These curves assume the low redshift luminosity-mass 
(specifically Lx(0A-2A keV) vs. Af2oo) relation from 
|Reiprich fc B5hringer| ( |2002| ) and, for simplicity, no 
redshift evolution. This scaling relation is also used 
to derive an X-ray mass estimate from the soft X-ray 
luminosity. Numerical values for these quantities are 
presented in the last two columns of Table [T] 
4. SUMMARY 

We have established that the central ellipticals associ- 
ated with the four SZE clusters have luminosities consis- 
tent with those of BCGs from the SDSS. We have also 
determined cluster masses using two different mass es- 
timators based on the luminosity of the optical galaxies 
and the X-ray-emitting gas. These mass tracers are rea- 
sonably consistent within the uncertainties given by the 
scatter in the X-ray luminosity-mass relation (~0.2 in 
the log) and the estimated factor of 2 uncertainty in the 
mas s derived from the opt ical luminosity of the galax- 
ies ( Menanteau et al~]|2008 |. Based on this evidence, all 
four clusters are fairly massive (M ^ 5 x 10 14 M Q ) and 
therefore well above the nominal detection limit of both 



flux end) that corresponds approximately to the Bright 

2 ftp:/ /ftp. xray.mpe.mpg.de/rosat/archive/ 



SPT and ACT (e.g., |Sehgal et al.||2007| ). Although this 
is an important step in identifying the optical and X- 
ray counterparts to the SZ clusters, it cannot be consid- 
ered definitive. In addition the RASS count rates could 
well be contaminated by X-ray emission from an AGN 
or other unrelated source. Thus the conservative reader 
could conclude that the X-ray luminosities and derived 
masses for the three higher redshift SPT sources repre- 
sent upper limits to their true values. Approved Chandra 
observations will within the next year resolve any doubts 
about the X-ray emission of these clusters. Deeper opti- 
cal imaging and spectroscopy should also be pursued in 
order to better constrain the masses of these clusters. 
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Note. — Physical properties of SZE selected clusters in the SCS regions. Redshifts represent the photometric redshift from the bright 
elliptical in the center of the cluster with ±ltr limits. The cluster position is based on the BCG. 
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